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The complex conductivity at the (Verwey) metal-insulator transition in Fe3C>4 has been inves- 
tigated at THz and infrared frequencies. In the insulating state, both the dynamic conductivity 
and the dielectric constant reveal a power-law frequency dependence, the characteristic feature of 
hopping conduction of localized charge carriers. The hopping process is limited to low frequencies 
only, and a cutoff frequency v\ ~ 8 meV must be introduced for a self-consistent description. On 
heating through the Verwey transition the low-frequency dielectric constant abruptly decreases and 
becomes negative. Together with the conductivity spectra this indicates a formation of a narrow 
Drude-peak with a characteristic scattering rate of about 5meV containing only a small fraction 
of the available charge carriers. The spectra can be explained assuming the transformation of the 
spectral weight from the hopping process to the free-carrier conductivity. These results support 
an interpretation of Verwey transition in magnetite as an insulator-semiconductor transition with 
structure-induced changes in activation energy. 
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I. INTRODUCTION 

Magnetite (Fe304) is probably the oldest known mag- 
netic material, which can be found in natural form. This 
material is of considerable importance because of var- 
ious applications in magnetic recording and for high- 
frequency electronic devices. Physical properties of mag- 
netite have attracted much attention after the discov- 
ery of the first-order metal-insulator transition at Ty ~ 
120 K. On cooling through the transition temperature the 
dc-conductivity drops by two orders of magnitude. A 
realistic model to explain the mechanism of this tran- 
sition has been suggested by Verwey 0, assuming a 
charge order-disorder transition with alternating valence 
(Fe 2+ /Fe 3+ ) of the octahedrally-coordinated iron ions. 
The metal-to-insulator transition in magnetite has been 
termed Verwey transition since then. In spite of a large 
number of experimental and theoretical efforts, the mech- 
anism governing the conduction and magnetic properties 
in magnetite is still under debate @,0|. According to the 
results of recent x-ray resonant scattering Q, even the 
charge-ordering below Ty has been questioned and the 
concept of an itinerant magnet was considered instead. 
Also the importance of the orbital degrees of freedom 
have recently been highlighted p|- 

Magnetite crystallizes in a cubic high-temperature 
structure and exhibits a monoclinic distortion below the 
Verwey transition. In spite of enormous progress in re- 
solving the low-temperature monoclinic structure of mag- 
netite 0, @ > full details could not be completely re- 
solved so far. The monoclinic phase is insulating and 
is characterized by a gap of the order of 100 meV in 
the excitation spectrum, which is seen in photoemission 
data [jj 0, E| , in thermoelectric properties , and by 



optical spectroscopy 0, The charge transport in 

ma gnet ite is usually explained within a polaronic picture 
However, there is no general agreement 
about the energy of the polaronic absorption in the con- 
ductivity spectra [la. X\M Ilia. Hq|. Recent photoemission 
results obtained in the same magnetite crystals as 
in the present work, have been self-consistently explained 
using a small-polaron model. A similar concept has been 
applied recently to explain the dynamic conductivity at 
the metal-insulator transition in quasi-one- dimensional 
/3-Na . 33 V 2 O 5 E3- 

At low-temperatures the conduction in magnetite takes 
place via hopping between localized states, which agrees 
well with both activated dc-resistivity 0, H^. l2fj| and 
with an observed characteristic power-law a oc v s pl| 
of the ac-conductivity. This universal power-law has 
been observed in various materials with internal disorder 
and is a clear fingerprint of hopping of charge carriers 
[23L I24I I2H . The situation is getting more complicated 
in the metallic state. In spite of two orders-of-magnitude 
increase in resistance above Ty, the absolute value of 
the dc-conductivity (a dc (T > Ty) ~ SOft^cm -1 ) is still 
much smaller than the Ioffe-Regel 26] minimum metal- 
lic conductivity a rn in ~ e 2 /3fia ~ 3000 51 _1 cm _1 , even 
including the Mott correction |2]| in case of disorder 
o'min ~ 0.03e 2 /?ia ~ 300f2 -1 cm -1 (here a ~ 3 A is the 
interatomic distance). The high-temperature conductiv- 
ity of magnetite is thus typical for a semiconductor with a 
high mobility of thermally-excited charge carriers rather 
than for a metal. Reviews of the recent and past devel- 
opments in magnetite can be found in 0, 0, IH, |2^, |3(j ■ 

In this paper we investigate Terahertz and infrared 
conductivity of magnetite on both sides of the metal-to- 
insulator transition. The unexpected result of these ex- 
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periments is the abrupt change of the dielectric constant 
at the transition temperature yielding even negative val- 
ues above Ty. Together with the conductivity spectra 
these results indicate the formation of a narrow band of 
quasi- free carriers which contains only a small part of the 
total spectral weight. 

II. EXPERIMENTAL DETAILS 

Synthetic single crystals of magnetite were prepared 
from a— Fe2 03 (spec, pure) using a floating-zone method 
with radiation heating j^H . Two single crystals were used 
for the measurements, which revealed slightly different 
transition temperatures, Ty ~ 123 ± 1 K and 116 ± 1 K, 
respectively. The transition temperature in magnetite 
can be taken as a criterium for impurity concentration 
and oxygen stoichiometry 2, 

which, therefore, differ 
significantly for the two crystals investigated. Both sam- 
ples revealed qualitatively similar electrodynamic proper- 
ties revealing only differences in absolute values. There- 
fore, in the following mainly the results on a sample with 
Ty ~ 123 K will be shown. Plane-parallel samples of dif- 
ferent thickness between 0.5 mm and 0.05 mm have been 
prepared from the original boules. To ensure mechan- 
ical stability, the thinnest samples were glued onto a 
MgO substrate of ~ 0.5 mm thickness. The availability 
of samples with different thicknesses is essential for trans- 
mittance experiments due to the strongly temperature- 
dependent absorption in magnetite. 

The dynamic experiments for frequencies 4 cm -1 < 
v < 40 cm -1 were carried out in a Mach-Zehndcr inter- 
ferometer arrangement |33J which allows both, the mea- 
surements of the transmittance and the phase shift of a 
plane-parallel sample. The complex transmission coeffi- 
cient has been analyzed using the Fresnel optical formu- 
las for transmittance T = \t\ 2 of a plane-parallel sample 

t = VTe^ = , (1) 

where r = and t\ — exp(— 2iTiy/e d/X) . 

V? + 1 

Here r is the reflection amplitude of a thick sample, t\ is 
the "pure" transmission amplitude, e is the complex di- 
electric permittivity of the sample, d is the sample thick- 
ness, and A is the radiation wavelength. It has been 
assumed that the magnetic permeability [i = 1 in the 
frequency range of our experiments. The transmittance 
for a two-layer system can be obtained in a similar way 
[34l Isflj . Using Eq. QJ the absolute values of the com- 
plex conductivity a* = u\ + 102 and dielectric constant 
e* = £\ +ie2 = cr* /ieoLu can be determined directly from 
the measured spectra. Here Sq and u = I-kv are the per- 
mittivity of vacuum and the angular frequency, respec- 
tively. In general, the phase shift <Pt in the transmittance 
experiment contains a possible ambiguity of 2-7T • n, where 
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FIG. 1: Left panels: temperature-dependent reflectance (top 
left) and reflectance phase-shift (bottom left) of magnetite at 
v = 12 cm -1 in the vicinity of the Verwey transition. Right 
panels: transmittance (top) and phase shift (bottom) of a 
0.05 mm thick plane-parallel sample glued on the 0.5 mm thick 
MgO substrate. The four measured quantities represent two 
data-sets to determine a\ and e\ = — 02 1 e^. 



n is an integer. To exclude this source of errors, the re- 
flectance of thick sample has been measured in addition 
to transmittance experiments. A similar Mach-Zehnder 
arrangement allowed to obtain the amplitude and the 
phase shift of the reflectance at the Verwey transition. 
In total, a set of four measured quantities was available 
to calculate the real and imaginary part of the dynamic 
conductivity in the THz frequency range. 

In the infrared and visible frequency ranges 30 cm -1 < 
v < 21000 cm -1 the complex conductivity has been ob- 
tained via the Kramers-Kronig analysis of the reflectiv- 
ity of the thick sample. The reflectivity experiment were 
performed using Brukcr IFS-113v and IFS 66v/S Fourier- 
transform spectrometers. Different sources, beam split- 
ters and optical windows allowed to cover the complete 
frequency range. In addition, the reflectance for the fre- 
quency range 4 < v < 40 cm" 1 has been calculated using 
the complex conductivity data, obtained by the trans- 
mittance technique, described above. The combination 
of the results from two experimental techniques |36j sub- 
stantially expands the low-frequency limit of the available 
spectrum and the quality of the subsequent Kramers- 
Kronig transformation. 



III. RESULTS AND DISCUSSION 

Figure ^ shows the complex reflectance and trans- 
mittance of magnetite in the vicinity of Verwey transi- 
tion. All four quantities reveal strong changes at the 
transition temperature. The most dramatic changes 
are observed in the transmittance, that decreases by 
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FIG. 2: Temperature dependence of the dielectric constant 
(ei, upper panel) and conductivity (<7i, lower panel) in mag- 
netite at v — 12 cm -1 . The data have been calculated from 
the complex reflectance (circles) of a thick sample and from 
the complex transmittance of a 0.05 mm thin plate (trian- 
gles). Right inset - dc-conductivity of a thick sample in full 
temperature range, left inset - temperature region close to the 
Verwey transition. 



three orders of magnitude. Both transmittance and re- 
flectance show a hysteretic behavior, characteristic for 
a first-order phase transition. However, the experimen- 
tally observed hysteresis extends over different temper- 
ature regimes: the hysteresis spans AT ~ 2 K in the 
reflectance experiments, and is significantly less than 
0.5 K in transmittance experiments. We suggest that 
this distinct difference probably results from the sample 
geometries, namely utilizing a thick sample (~ 5 mm) 
in the reflectance and thin sample (~ 0.05 mm) in the 
transmittance experiments. Except for the hysteresis re- 
gion, all measurements revealed results coinciding within 
the experimental accuracy (Fig-EJ) . The effective over- 
determination of the physical quantities, £\ and <7i, in 
these experiments was necessary because of the highly 
unusual temperature and frequency dependence of the 
electrodynamic properties of magnetite (as discussed be- 
low). 

Figure [21 shows the temperature dependence of the 
conductivity (o\) and the dielectric constant (e\) of 
magnetite at v = 12 cm -1 . The data have been ob- 
tained either from the complex transmittance (t = 



VT exp(i(f>r) , triangles) or from the complex reflectance 
(r = \/Rexp(i(j)ii), circles). The conductivity reveals a 
sharp increase by a factor of three at the transition tem- 
perature. However, this value is substantially lower than 
the increase by two orders of magnitude, observed in the 
dc-conductivity (insets of Fig.[21> ■ This is a direct con- 
sequence of the power-law frequency dependence of the 
ac-conductivity yielding orders-of-magnitude increase of 
the conductivity with increasing frequency. 

The behavior of the dielectric constant at the Verwey 
transition is highly unusual. Both transmittance and re- 
flectance experiments lead to a strong decrease of the 
dielectric constant at the metal-to-insulator transition, 
becoming even negative above Ty. In the classical sys- 
tems revealing a metal-to-insulator transition (e.g. doped 
silicon) it has been observed that the dielectric constant 
diverges approaching the transition temperature from the 
insulating side |27l l37| , which has been termed dielectric 
catastrophe. In agreement with this picture and using the 
same technique as presented here, a strong increase of the 
dielectric constant at the metal-insulator transition has 
been observed in Sr-doped LaMnC>3 |38| and recently in a 
magnetic field-induced transition in Pr(Ca:Sr)MnC>3 [39|. 
We believe that the first-order character of the phase 
transition in magnetite interrupts the increase of the di- 
electric constant, observed at T < Ty (de/dT ~ 1 K _1 at 
low frequencies). Instead, the dielectric constant of mag- 
netite abruptly jumps to negative values at the "metallic" 
side of the transition. 

Figure [21 shows the frequency dependence of the con- 
ductivity and dielectric constant in magnetite. Symbols 
at low frequencies represent the results of the trans- 
mittance experiments. Solid lines above 30 cm" 1 have 
been obtained via the Kramers-Kronig analysis of the re- 
flectance. The lower frame of Fig. [3] shows the frequency 
dependence of <j\ above and below the Verwey transition. 
The low-frequency conductivity in the insulating state is 
dominated by a power-law in frequency (&\ oc v s with 
s 1.3) with a temperature-dependent amplitude and a 
weakly temperature-dependent frequency exponent. As 
mentioned in the Introduction, a power-law behavior of 
the conductivity is a characteristic feature of hopping 
conduction between localized states [23, 0j HH • Previ- 
ously, the frequency exponent s ~ 0.7 has been observed 
in magnetite at low temperatures and at kHz-frequencies 
|2l| , a value typical for the ac-conductivity in the audio- 
frequency range. Substantially higher values of the ex- 
ponent are expected at higher frequencies: approaching 
the phonon-assisted-hopping regime, <7i oc v 2 has been 
predicted for low temperatures 40] . The transition to 
higher frequency-exponents in the microwave frequency 
range has been discussed recently for doped semiconduc- 
tors |4l| . The power-law term is g enerally referred to 
as universal dielectric response .22| and has been a re- 
cent matter of discussion concerning a possible universal 
superlinear power law with s > 1 at high frequencies |42| . 

The frequency dependence of the dielectric constant 
in magnetite is shown in the upper panel of Fig. [3 
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FIG. 3: (color) Frequency dependence of the conductivity (<7i, 
lower panel) and the dielectric constant (ei, upper panel) 
of magnetite above and below TV — 123 K. Symbols below 
40 cm -1 and solid lines above 30 cm -1 represent the experi- 
mental data, dashed lines - model calculations. Experimental 
data have been obtained from transmittance and phase-shift 
of a thin FesCU plate below 40 cm -1 , and via the Kramers- 
Kronig analysis of the reflectance of the thick sample above 
30 cm - . The main contribution to the model calculations is 
given by the hopping term (a* cx v s , Eq.|2J below TV and a 
Drude-term (Eq.0J above Ty. 



The behavior of the dielectric constant in the insulat- 
ing state agrees qualitatively with the power-law of <7i, 
i.e £\ increases with decreasing frequency. For s > 
1 the low-frequency behavior can be approximated by 
ei(0) - £i(» oc V s ' 1 . 

Dashed lines in Fig.[3]for the insulating state (T < Ty) 
were obtained by fitting the experimental data to the 
expression 



<Ti(v) =a dc + A 



1 + {v/vi) A 



(2) 



Here a<ic is the dc-conductivity and A is the amplitude of 
the hopping process. Compared to the conventional form 
of this response [H El Si IU E2, we introduced an 
additional high-frequency cutoff [1 + [y fv\f\ ■ The cutoff 
frequency v\ only weakly influences the conductivity in 
the low-frequency range but prevents the divergence of 
the spectral weight of o\ cx v s at high frequencies. 

The analytical expression for the dielectric constant 
has been obtained by applying the Kramers-Kronig 
transformation |4i| to Eq. which gives the following 



TABLE I: Parameters of the THz-frequency excitations in 
magnetite close to the Verwey transition which correspond 
to Eqs. 1213 14H . The changes in conductivity can be well de- 
scribed assuming the transformation of the hopping term 
Eq. into the Drude-term Eq. JIJ with a high mobility of 
the quasi-free carriers. 



T (K) vx (cm" 1 ) 


s 


A(fi -1 cm _1 ) 


100 70 


1.5 


39 


120 60 


1.2 


44 


T (K) l/27rr(cm- 1 ) 




cr (ic (n~ 1 cin -1 ) 


130 40 
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analytical expression for the dielectric constant: 



cos(s7r/2) 



<Ti(z/) ~ Ode) , , /r) , 

£ i( v ) = — 7, i tan(s7r/2) 



sm[0 + l)-] + [ — J cos[(« + l)~] 



(3) 



This expression is valid both for s < 1 and for s > 1. 

To take into account the phononic and high-frequency 
electronic contributions, the sum of 6 Lorentzians has 
been added to expressions Eq. (J2J) and Eq. @ . Five 
Lorentzians at frequencies of 106, 152, 247, 361, and 
574 cm -1 represent phonon contributions. Additional 
over-damped oscillator at 4270 cm -1 (0.53 eV) approxi- 
mates the first excitation of the electronic origin. The 
given characteristic frequencies agrees well with the pub- 
lished data [Til. FLU Il5l | . In order to simplify the analysis 
of the THz conductivity and obtain the essential physics 
of the problem, the contributions from the Lorentzians 
were fixed for all three temperatures in Fig-EI The pa- 
rameters obtained from simultaneous fitting of the con- 
ductivity and dielectric constant are given in TableQ] 

A remarkable feature of the parameters for the low- 
temperature insulating phase in Tabled is the low value 
of the cutoff frequency v\. Indeed, a qualitative exami- 
nation of the conductivity frame of Fig.[3]reveals that the 
hopping process o\ oc v s and the low-frequency edge of 
the electronic excitation [y > 1000 cm -1 ) are well sepa- 
rated in energy. This indicates that the cutoff frequency 
of the hopping must be in the range ~ 100 cm , in agree- 
ment with the fits. This result is rather surprising and 
suggests that the low-temperature conductivity is gov- 
erned by phonons. 

In contrast to the gradual variation of £i(o>, T) and 
oi(u),T) at low temperatures, dramatic changes are ob- 
served at the Verwey transition. At low frequencies 
the conductivity increases by more than one order-of- 
magnitude. This increase is substantially smaller than 
a change of two orders-of-magnitude in dc-conductivity 
(see insets of Fig.[2Jl. This difference is the result of 
a power-law frequency dependence of the conductivity 
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in the insulating state. As a consequence, for frequen- 
cies close to 35 cm -1 we observe almost no conductivity 
change at Ty • This observation has been verified in a sep- 
arate temperature-dependent experiment at v — 35 cm -1 
(not shown). 

As discussed above, the decrease of the dielectric con- 
stant at the transition to the conducting state is not ex- 
pected within the scope of a conventional metal-insulator 
transition scenario. In order to find an explanation to 
this experimentally-observed behavior, we mention that 
a negative dielectric constant naturally follows by assum- 
ing a Drude- response of quasi-free carriers, 



'DR. 



0"dc/(l ~ ilnvr) 



(4) 



where <Jdc and 1/2ttt are dc-conductivity and scattering 
rate, respectively. In addition to the negative changes in 
£i, the conductivity in the metallic state shows a slight 
downward curvature in the frequency dependence (lower 
panel of Fig.[3|). This provides a further argument to 
include a Drude-term to the model. Indeed, the simplest 
way to reproduce the experimental changes between 120 
and 130 K is just to replace the hopping term fEa. l2l3|) 
by the Drude term. For the second sample of magnetite 
with Ty — 116 K a qualitatively similar behavior has 
been observed and the absolute values of the parameters 
agreed within ~ 20%. We believe that this variation 
partly corresponds to the difference in the conductivities 
of the samples. 

The low value of the scattering rate, as observed in 
the present experiments, corresponds to a high carriers 
mobility of ji = er/m ~ 10 3 cm 2 /ys (assuming m be- 
ing the free electron mass) and correlates well with the 
high quality of our magnetite single crystals. Due to the 
low characteristic frequency of the coherent (itinerant) 
process, its spectral weight 



S = — J ai. DR (uj)duj = o-dc/r 
o 



(5) 



is small compared to the full number of carriers. Tak- 
ing as an estimate of the total spectral weight the con- 
centration of Fe 2+ ions j2(|, only about 2 • 10~ 4 of the 
theoretically-available carriers participate in the coher- 
ent process. We recall further, that within the rough es- 
timate the spectral weight of thermally-activated quasi- 
free carriers in a semiconductor can be written as 



S = 



l eff 



■ exp(- 



A, 



(6) 



where n, e and m e ff are concentration, charge and ef- 
fective mass of the charge carriers, respectively. In this 
expression the effective concentration of charge carriers is 
temperature-activated and is governed by an energy gap 
A m . The estimate of the effective mass m e // — 10 2 m 



has been obtained from the infrared conductivity |l4| and 
from the small-polaron analysis of photoemission of the 
same sample [Tl|. With this estimate for m e // the ob- 
served small spectral weight of the Drude-process corre- 
sponds to A m ~ 40 meV, which correlates with a weakly- 
activated dc-conductivity just above Ty (right inset in 
Fig.EJ. Within the same picture, the Verwey transition 
in magnetite corresponds just to an increase of the char- 
acteristic activation energy at the insulating side of the 
transition. The drop of dc-conductivity at the transi- 
tion can be ascribed to a gap change from A m ~ 40 meV 
(conducting state) to A, ~ 90meV (insulating state). 
This new value of the energy gap agrees with the ac- 
tivation energy of dc-conductivity 0, 0, > thermo- 
electric power |12j . and with the gap values obtained in 
photoemission experiments |9t ll0llll| . The presented de- 
scription implies the interpretation of Verwey transition 
as an insulator-semiconductor transition with a change 
in activation energy induced by crystal structure or by 
electronic configuration. 



IV. CONCLUSIONS 

Terahertz and infrared conductivity in magnetite has 
been obtained on both sides of Verwey metal-to-insulator 
transition. Two orders-of-magnitude change in dc- 
conductivity halves at millimeter frequencies and even 
disappears for v > 30 cm . In the far-infrared range 
the effect of the transition remains solely in changing the 
dielectric constant (ei) of magnetite. The change of e\ 
during the transition into the conducting state is nega- 
tive, which is in surprising contrast to an expected (pos- 
itive) divergence of the dielectric constant at a metal- 
insulator transition. Within the simple model-analysis 
the conductivity mechanism switches between hopping 
of localized carries below Ty and itinerant motion above 
Ty. These results evidence the formation of the coher- 
ent state on the metallic side of the Verwey transition 
with high mobility of the charge carriers. Together with 
the recent photoemission data on the same samples the 
interpretation of the Verwey transition as an insulator- 
semiconductor transition is suggested. 
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